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The diffusive behavior of a harmonic oscillator driven by a Mittag-Leffler noise is studied. Using 
Laplace analysis we derive exact expressions for the relaxation functions of the particle in terms of 
generalized Mittag-Leffler functions and its derivatives from a generalized Langevin equation. Our 
results show that the oscillator displays an anomalous diffusive behavior. In the strictly asymptotic 
limit, the dynamics of the harmonic oscillator corresponds to an oscillator driven by a noise with a 
pure power-law autocorrelation function. However, at short and intermediate times the dynamics 
has qualitative difference due to the presence of the characteristic time of the noise. 
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1. INTRODUCTION 

The study of anomalous diffusion in complex or disor- 
dered media has achieved a substantial progress during 
the last years [H, 0, H, 0, [E B ■ Anomalous diffusion in 
physical and biological systems can be formulated in the 
framework of the gene ralized Langevin equation (GLE) 
[H, 0, H S, ES mill El EI O- If one considers the 
dynamics of a harmonic oscillator with frequency to un- 
der the influence of a random force modeled as Gaus- 
sian colored noise, the corresponding GLE is written as 



X(t) + / dt' 1 {t-t')X{t')+uj 2 X = £{t) ) 



(1) 



where X(t) represents the position of a particle of mass 
772 = 1 at time t, and "/(t) is the frictional memory kernel. 
The random force is zero-centered and stationary 
Gaussian that obeys the fluctuation-dissipation theorem 



mm} = c(\t-t'\) = k B T 1 (\t-ti\), 



(2) 



where ks is the Boltzmann constant, and T is the abso- 
lute temperature of the environment. 

It is now well established that the physical origin of 
anomalous diffusion is related to the long-time tail cor- 
relations [l|, [2|. Therefore, in order to model anomalous 
diffusion process, pure power-law correlation functions 
are usually employed [J B El El El E! , 

which may be 

written as 



C(t) = C x 



\t\ 



r(i-A)' 



(3) 
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where T(z) is the Gamma function, and C\ is a pro- 
portionality coefficient dependent on the exponent A but 
independent of time. The exponent A can be taken as 
0<A<lorl<A<2, which is determined by the dy- 
namical mechanism of the physical process considered. 

Vinales and Desposito have introduced a noise whose 
correlation is proportional to a Mittag-Leffler function 
[2Cj |. This correlation behaves as a power-law for large 
times, but is non-singular at the origin due to the inclu- 
sion of a characteristic time. 

The aim of this work is to investigate the effects of the 
Mittag-Leffler noise on the behavior of a harmonically 
bounded particle governed by the GLE JT]). This paper 
is organized as follows. In Section 2 we discuss some 
characteristics of the Mittag-Leffler noise. In Section 3, 
we show the formal expressions for the relaxation func- 
tions that govern the dynamics of the particle in the case 
of an arbitrary noise correlation function. Analytical so- 
lutions of the GLE for a harmonically bounded particle 
driven by a Mittag-Leffler noise are obtained in Section 
4. The Section 5 is devoted to the analysis of tempo- 
ral behavior of the relaxation functions, and is compared 
with the results in the case of a pure power-law noise cor- 
relation function. Finally, the conclusions are presented 
in Section 6. 



2. MITTAG-LEFFLER NOISE 

It is well known that if the correlation function ^ is a 
Dirac delta function the stochastic process is Markovian 
and its dynamics can be directly obtained [2l|. However, 
in a complex or viscoelastic environment one must take 
into account the memory effects through a long-time tail 
noise to describe the effect of the environment on the par- 
ticle. The non-Markovian dynamics is involved in these 
physical processes. 

Recently, Vinales and Desposito introduced a Mittag- 
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Leffler noise given by [2C 
C(t) = % 



where 



E x (-(\t\/r) x ), 



(4) 



where r acts as a characteristic memory time and < 
A < 2. The E a (y) function denotes the Mittag-LefHer 
function [22j defined through the following series 



E a (y) 



i) 



a > 0. 



(5) 



Using the asymptotic behaviors of the Mittag-LefHer 
function [23[ one can easily deduce that, for A ^ 1, the 
correlation function (0| behaves as a stretched exponen- 
tial for short times and as an inverse power law in the 
long time regime 0, [24| • 

Setting A = 1, the correlation function ((4]) reduces to 
an exponential form 



T 



t|/7 



(6) 



which describes a standard Ornstein-Uhlenbeck process 
(2l| . Moreover, in the limit r — > and from the 
limit representation of the Dirac delta [25| we get that 
C(t) = 2 C\ 5(t) which corresponds to a white noise, non- 
retarded friction and standard Brownian motion [211 ] . 

On the other hand, for A ^ 1 the limit r — ► of the 
proposed correlation function (]4]) reproduce the power- 
law correlation function This behavior is obtained 
introducing in expression (|4]) the asymptotic behavior at 
large y of the Mittag-LefHer function (23| 



E a (-y) ~ [yT(l- a )Y 



y>0. 



(7) 



It is worth pointing out that the Mittag-LefHer cor- 
relation function Q is a well defined and no-singular 
function. From its value at t — is given by 
C(0) = C\/t x , while for the power-law correlation ^ 
C(0) diverges. Then, the introduction of the characteris- 
tic time r enables to avoid the singularity of the power- 
law at the origin. Considering that the Mittag-LefHer 
function is the natural generalization of the exponential 
function [12] , we can also consider the Mittag-LefHer cor- 
relation function as a generalization of the power-law cor- 
relation, and similarly, the colored noise ^ is considered 
as a generalization of the white noise. 



3. SOLUTIONS OF THE GENERALIZED 
LANGEVIN EQUATION 

In what follows we consider the Langevin equation (fT]) 
with the deterministic initial conditions xq = X(0) and 
t>o = -^(0). By means of the Laplace transformation to 
Eq. (TT]) one can easily obtain a formal expression for the 
displacement X(t) and the velocity V(t) — X(t). The 
displacement X{t) satisfies that 



X(t) = (X{t)) + / dt'G(t - t')£{t') , 



(8) 



(X(t)) =v G(t)+xo{l-uj 2 I(t)) 



(9) 



is the position mean value. The relaxation function G(t) 
is the Laplace inversion of 



G(s) = 



1 



s 2 + 7( s ) s + ^ 



(10) 



where j(s) is the Laplace transform of the damping ker- 
nel and 



I(t)= / dt'G{t') 



(11) 



On the other hand, the velocity V(t) satisfies that 



V(t) = (V(t)) + / dt'g(t - t'W) , (12) 
Jo 



where 



(V(t)) =v g(t)-Lo 2 x G(t) 



(13) 



is the velocity mean value and the relaxation function 
g(t) is the derivative of G(t), i.e. 



g(t) = G'(t). 



(14) 



Explicit expressions of the variances can be obtained 
from Eqs.© and (|12[) . Taking into account the symme- 
try property of the correlation function and Eq.([2]), yields 

f3a xx (t) = 2I{t)-G 2 (t)-u; 2 I 2 {t), (15) 
(3a vv (t) = l-g 2 {t)-u 2 G 2 {t) 1 (16) 
pa xv (t) = G{t){l-g{t)~io 2 I{t)} , (17) 

where /3 = l/k B T. 

From an experimental point of view, information about 
the observed difiusive behavior is extracted from the 
mean square displacement p{t). In the long time mea- 
surement time, it is related to the relaxation function 
I(t) as [H] 

P (t l ) = lim ((X(t + tl) - X{t)f) = 2k B TI( TL ) , (18) 

t — >oo 

where tl is the so-called time lag. Alternative informa- 
tion about the dynamics can be extracted from the nor- 
malized velocity autocorrelation function Cy (t), which is 
related to the relaxation function g(t) as [16, 26] 

Then, the knowledge of the relaxation functions I(t), 
G(t) and g(t) allows us to describe the difiusive behavior 
of the oscillator. In the next section we will give explicit 
expressions for the relaxation functions in the case of a 
Mittag-LefHer noise (01 assuming that A ^ 1. 
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4. ANALYTICAL RELAXATION FUNCTIONS 
FOR A MITTAG-LEFFLER NOISE 

From relation |(2J), the memory kernel j(t) correspond- 
ing to the Mittag-Leffler noise (U|) can be written as 



7 (t) = ^E x (-(\t\/ry) 



(20) 



where 7a = C\/ksT . Taking into account that the 
Laplace transform of the memory kernel reads [23[ 



7(s) 



7as 



A-l 



1 + S A T 



A—A 



(21) 



the relaxation function I(t) can be written as the Laplace 
inversion of 



j( a )= <M =I ( s )+h(s) 



where 



Io(s) 



T A S 2+A _|_ s 2 _|_ <-y AS A _|_ ^2 



h(s) = T A s A f ( S ), 
and 7a is defined as 



7a = 7a + u 2 t x 



(22) 

(23) 
(24) 

(25) 



Following the approach given in Ref. [ 2Jj we get 



n\ ^-^ ml 

m=0 



x« 2 4? + t + A + ( 2 -A )m (-(*/r) A ), (26) 



where To(^) an d -M*) are given by ((26)) and (|27f respec- 
tively. 

The relaxation functions G(t) and ffft) can be calcu- 
lated using |TTJ, lHU) and the relation (27| 



Then, we get 



where 



G (t) 



and 



where 



9o{t) 



G(t) = G (t) + G 1 (t) 



(31) 



(32) 



n=0 



A 6 ^A,2+2n+A+(2-A)m 



m=0 
H*A0 



(33) 



X ^£ + + t + (2-A, m (-(^) A ) 



fl(*)=S0(«)+Si(*) 



(34) 



(35) 



Tj l Z—^ 77-7 I 



^ n! 

n— m— 

X ^A?i+2«+A+(2-A)m( _ ( < / T ) A ) : 



(36) 



^e^V^eW- 



n— m— 

.2 ^(n+m) 



(27) 



where E a .p(y) is the generalized Mittag-Leffler function 
[23| defined by the series expansion 



£a,0(tf) = E 



y J 



3=0 



T(aj + (3) 



, a>0, /3>0, (28) 



and E^ g(y) is the derivative of the Mittag-Leffler func- 



tion 



^) = i^) = E 



Then, from 



I(t)=I (t)+h(t). 



(29) 



(30) 



5i W 



OO 

E 

n=0 



00 -2>41 

y \ rA 



x^T + t +(2 -A )m (-(Vr) A ). 



(37) 



It is worth mentioning that expressions (|30"|) . (|32| and 
(|35[) fully determine the temporal evolution of the mean 
values © and (fl"3"|) , variances (fl"5)) to (TIT)) , mean square 
displacement (| L8[) and velocity autocorrelation function 

Notice that in the limit w — > only survive the terms 
with n = in equations (26|) and l[27]). Then, Eq. ([25]) 
reduces to 7a = 7a, and the expression of the relaxation 
function I(t) for the free particle case is recovered. 

On the other hand, in the limit r — > the function 
2i(i) vanishes and the behavior of Io(t) can be achieved 
introducing the asymptotic behaviors of the generalized 
Mittag-Leffler function [27| 



E a A~y) 



1 



yT(P-a)- 



y >0, 



(38) 
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and its derivative 



and 



E ( a k U-y) 



1 



J Q ,/3V vr- yk +i T(P-a) 
in Eq. (|26|) . Then, after some algebra we obtain 



(39) 



I(t) = lim J„(t) = ^ 

n=0 



X* 2 ^,3 + A n (-7At 2 - A ) 



(40) 



-A oo 

E 

n=0 



(-1)" 



2,2\n 



x i 2 K (n) 



(-7At 2 - A ) 



2-A,3-A+An 



(47) 



Let us analyze the behaviors of the relaxation func- 
tions I(t), G(t) and git) for 7a £ 2 ~ A ^> 1. Introducing 
the approximation ([39]) in (j46|) and (f47|) . after some cal- 
culations and using §5§ one gets 



where we have used that 7\ — > 7a for r — > 0, according to 
([25| . The expression in series given in Eq. (|4"0|) coincides 
with the expression for the relaxation integral function 
I{t) corresponding to a pure power-law correlation func- 
tion, previously obtained in Ref . (l6j . 

Likewise, one can verify that in the limit r — > the 
relaxation functions G{t) and g(t) are also the same to 
that in the case of a pure power-law correlation function, 
given by 



G(t) 
9(t) 



E 

00 

E 

n=0 



(_, ,2+2\ n 



, ,2 +2\ n 

1 ) P (n) / 
I ^2-A,l+AnV 



-7A* 



2-A\ 



(42) 



1(f) w-^r U-veA-v— t 
I V 7a 

Then, from ([11]) and (fT?]) we get 



v d „ / w 



and 



" ^ r, / ^ A 

d 2 f y 7A 



where we introduced the dimensionless factor 
v=?±, < v < 1 . 

7A 



(48) 



(49) 



(50) 



(51) 



5. TEMPORAL BEHAVIOR OF THE 
RELAXATION FUNCTIONS 

The analytical expressions (|5D|) . and ([33)1 are the 
main result of this work. In the following we will analyze 
the time behavior of the relaxation functions for different 
regimes. 

The short-time behavior (t <C r) of the relaxation func- 
tions can be obtained using the series expansions (|28[) and 
(1291). Then 



I(t) 
G(t) 
5(0 



2 

* - 
1 - 



T A 



*_ 

24' 



/7a 



2a 

-r A 



f 2 
2 ' 



(43) 
(44) 
(45) 



which are the expected for a harmonic oscillator driven 
by a noise with a finite correlation at the origin [jjl. Hoi. |28| . 

Now we get an expression of the function lit) for times 
bigger than the characteristic time r of the noise, i.e. 
Or. For this purpose we introduce the approximation 
(|3"9"| in Eqs. (J2U) and (f2"T)l . After some algebra we get 



n=0 



( W V) r 



x^^-WaJ^a^), (46) 



The relaxation functions (|4"8)) to ([50)1 have the same 
functional form to the results obtained in the pure power- 
law case [lH but with the presence of the scale factor 
v. In the limit r — > is 1^ = 1 and one recovers the 
expressions corresponding to a pure power-law noise [lj| . 

It is worth pointing out that these expressions are the 
same to those that can be obtained directly discarding 
the inertial term s 2 in fT0"|). Then, Eqs. |[2SJ) to (5DJ) 
represents the solutions in the high friction limit. 

The strictly asymptotic behavior of the relaxation 
functions lit), G(t) and git) can be obtained introducing 
the asymptotic behavior of the Mittag-Leffler func- 
tion in Eqs. f4"8" l) to ([50 ] ). Then, for v ^ t x > 1 the 
relaxation functions can be written as 



1 7 A sin(ATr) T(A) 

U! 2 U! 4 7T t X 

7a sin(ATr) T(A + 1) 

CO 4 IT t A +! ' 

_7A sin(AT r) T(A + 2) 

LU 4 7T 



t \+2 



(52) 

(53) 
(54) 



lit) 
Git) 

g(t) 



As expected, the relaxation functions (|52|) to ([54)1 be- 
have as a power law in the long-time limit. These results 
are in agreement with those obtained in Refs. [16l. |29|| due 
to the fact that the Mittag-Lefler noise decays as a power- 
law for very large times. In the same way, substitution of 
these asymptotic expansions into Eqs. (TT5")) to (fTT)) give 
the long-time behavior of the variances of the pro cess, 
which again coincide with those obtained in Refs. |l6l.l29l|. 
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6. CONCLUSIONS 

In this work we have presented an analytically resol- 
uble model for the dynamics of a classical harmonic os- 
cillator in a complex environment, which is valid for all 
time range. We have shown that an anomalous diffusion 
process can be generated by a Mittag-Leffler noise de- 
riving exact expressions for the relaxation functions of 
the oscillator in terms of the generalized Mittag-Lcfflcr 
function and its derivatives. Moreover, in the appropri- 
ate limits the results for a harmonic oscillator driven by 
a power-law noise is recovered. However, differences in 
relation to the usually employed pure power-law noise ap- 
pear in the interval of short and intermediate times. For 
times shorter than the characteristic time of the noise 
the relaxation functions include a correction due to the 
presence of the characteristic time r. In the range of 
intermediate times, the relaxation functions have a simi- 
lar functional form to the previously obtained for a pure 
power-law noise [Tf3], but with the inclusion of a scaling 
dimensionless parameter. Finally, in the strictly asymp- 



totic limit we recover the anomalous behavior of an har- 
monically bounded particle driven by a power-law noise, 
which is in agreement with the previous results given in 
Refs.HEHI. 



Acknowledgments 

This work (Desposito and Vihales) was performed un- 
der Grant N° PICT 31980/05 from Agenda Nacional dc 
Promotion Cienti'fica y Tecnologica, and Grant N° X099 
from Universidad de Buenos Aires, Argentina. K.G. 
Wang is pleased to acknowledge partial financial sup- 
port received from the Materials World Network Pro- 
gram and Metallic Materials and Nanostructures Pro- 
gram of National Science Foundation, Washington, DC, 
under Grants DMR-0710484 and from Florida Solar En- 
ergy Center. K.G. Wang is grateful to Drs. Desposito 
and Vihales' warm hospitality and fruitful discussions 
when he visited there in 2008. 



K. G. Wang and C. W. Lung, Phys. Lett. A 151, 119 [16 
(1990); K. G. Wang, Phys. Rev. A 45, 833 (1992). 
R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000); R. [17 
Metzler and J. Klafter, J. Phys. A: Math. Gen. 37, R161 
(2004). [18 

A. Bunde, S. Havlin, J. Klafter, G. Graff, A. Shehter, [19 
Phys. Rev. Lett. 78, 3338 (1997). [20 
E. Barkai and R. J. Silbey, J. Phys.Chem. B 104, 3866 
(2000). [21 
E. Barkai, Phys. Rev. E 63, 046118 (2001). 

P. Dieterich, R. Klages, R. Preuss and A. Schwab, Proc. [22 
Natl. Acad. Sci. U.S.A. 105, 459 (2008). 

B. Yilmaz, S. Ayik, Y. Abe and D. Boilley, Phys. Rev. E [23 
77, 011121 (2008). 

J. M. Porra, K. G. Wang and J. Masoliver, Phys. Rev. E [24 
53, 5872 (1996). [25 
K. G. Wang and J. Masoliver, Physica A 231, 615(1996). 
K. G. Wang and M. Tokuyama, Physica A 265, 341 [26 
(1999). [27 
S. Chaudhury, D. Chatterjee and B. J. Cherayil, J. Chem. 
Phys. 129, 075104 (2008). [28 
S. Chaudhury, S. C. Kou, and B. J. Cherayil, J. Phys. 
Chem. B, 111, 2377-2384 (2007). [29 
S. C. Kou, X. S. Xie, Phys. Rev. Lett. 93, 180603 (2004). 
K. S. Fa, Phys. Rev. E 73, 061104 (2006). 
K. S. Fa, Eur. Phys. J. B 65, 265-270 (2008). 



A. D. Vifiales and M. A. Desposito, Phys. Rev. E 73 
016111 (2006). 

R. Zwanzig, Noneqmlibrium Statistical Mechanics (Ox- 
ford Univ. Press, New York, 2001). 

E. Lutz, Europhys. Lett. 54, 293 (2001). 

S. Burov and E. Barkai, Phys. Rev. E 78, 031112 (2008) 
A. D. Vifiales and M. A. Desposito, Phys. Rev. E 75 
042102 (2007). 

H. Risken, The Fokker-Plank Equation, (Springer- Verlag 
Berlin, 1989). 

A. Erdelyi et al, Higher Transcendental Functions. 
(Krieger, Malabar, 1981), Vol.3. 

F. Mainardi and R. Gorenflo, J. Comput. and Appl 
Mathematics 118, 283 (2000). 

R. Granek and J. Klafter, Europhys. Lett. 56, 15 (2001) 
C. Cohen- Tannoudji et al, Quantum Mechanics ,Vol. II 

(1977). ' 

M. A. Desposito and A. D. Vifiales. larXiv:0902.278"6V l. 

I. Podlubny, Fractional Differential Equations (Aca- 
demic Press , London, 1999). 

A. D. Vifiales, Ph. D. Thesis, Universidad de Buenos 
Aires, (2008). 

M. A. Desposito and A. D. Vifiales, Phys. Rev. E 77, 
031123 (2008). 



